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ABSTRACT 

In this study, we present a method for optimizing, at the 
preliminary design level, a supersonic turbine tor rocket 
propulsion system application. Single-, two- and three- 
stane turbines are considered with the number ol design 
variables increasing from b to 11 then to 15. in 
accordance with the number of stages. Due to its global 
nature and flexibility in handling different types ol 
information, the response surface methodology (RSM) is 
applied in the present study. A major goal of the present 
optimization effort is to balance the desire of maximizing 
aerodynamic performance and minimizing weight. To 
ascertain required predictive capability of the RSM. a 
two-level domain refinement approach has been adopted. 
The accuracy of the predicted optimal design points based 
on this strategy is shown to be satisfactory. Our 
investigation indicates that the efficiency rises quickly 
from single stage to 2 stages but that the increase is much 
less pronounced with 3 stages, A I -stage turbine performs 
poorlv under the engine balance boundary condition. A 
significant portion of fluid kinetic energy is lost at the 
turbine discharge of the 1 -stage design due to high stage 
pressure ratio and high-energy content, mostly hydrogen, 
of the working fluid. Regarding the optimization 
technique, issues related to the design of experiments 
<DOF) has also been investigated. It is demonstrated that 
the criteria for selecting the data base exhibit significant 
impact on the efficiency and effectiveness ot the 

construction of the response surface. 

NOMENCLATURE 

A, nil l:\it Blade Annulus Area, 
c : Chord 

D : Diameter 

h : Blade Height 

■ . jnii-m 1 2U00 h\ 1'uhlisiicil in the Amencaii institute ot AaoiutiU-.n .iiu 

Astronautics. Inc. with permission* 


fii : Mass flow 

P : Pressure 

R : : Proportion of Variation 

nns : Root Mean Square 

RPM : Angular Speed 

r : Standard Deviation 

se : Standard Error 

^r : Stage Reaction 

T n : Input Temperature 

: Pitch Speed 
\V : Weight 

a f : Work Fraction 

Apay : Payload Increment 

q Efficiency 

L INTRODUCTION 

Supersonic turbine technologies are being actively 
investigated in the rocket propulsion community. 
Optimizing a multistage turbine is a labor-intensive task 
and it is desirable to develop efficient and effective 
techniques to undertake this task. In general, two types of 
optimization are needed, namely, 

l Preliminary design* in which Mmplified models 
employing loss correlations gleaned from 
experimental database and one-dimensional 
gasdynamic and thermodynamic considerations, and 

2. Detailed shape design of the turbine blades , in which 
three dimensional computational fluid dynamics and 
detailed experimental information is employed. 

In this study, we present an approach based on the 
response surface methodology 1 (RSM) for optimizing, at 
the preliminary design level, a supersonic turbine aimed 
for Reusable Launching Vehicle (RLV) propulsion 
system application. Single-, two- and three-stage turbines 
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are considered with the number of design variables 
increasing in accordance with the number of stages. In the 
first step, the overall dimensions such as mean diameter 
and axial chords, the RPM and the number of stages are to 
be determined. In the second step, detail blade geometries 
will be optimized to achieve the best efficiency for a 
given overall gas path meridional geometry. This paper is 
limited to the preliminary optimization. 

In the past such optimization tasks might take weeks to 
perform due to the large number of parameters involved. 
The systematic application of RSM computationally 
coupled with an appropriate turbine analysis code will in 
the future allow designers to cut this cycle time down to a 
few hours. 

There are 2 types of design variables: 

1. Geometric input needed to layout the turbine 
meridional geometry. 

• mean diameter 

• last rotor annulus area 

• blade height ratio between the P vane and the 

last rotor blade - linear distribution of blade 
heights is assumed between the P vane and the 
last rotor blade 

• vane and blade axial chords 

2. Performance input used to calculate the turbine 
efficiency 

• RPM 

• number of stages 

• blade row reaction 

• work split (if more than 1 stage is investigated) 

The above parameters are essentially diffident to 
determine the geometry ,md the peri orma nee oi a 
preliminarv design. As indicated in tables 1-3. lor 
turbine with 1 stage, 2 stages or 3 stages the number ot 
parameters are 6, 11, 15, respectively. Constraints are also 
part of the optimization process. There are 2 structural 
constraints, the blade centrifugal stress and the disk stress. 
The blade centrifugal stress was constrained by a limit 
placed on the lumped inertia measure (the product of the 
blade exit annulus area and the RPM~ ). The disk stress 
was constrained by a limit placed on the pitch line velocity 
(the product of the RPM and the mean radius). 

For rocket engine applications, maximizing the vehicle 
payload for a given turbine operating condition is (he 
ultimate objective. Any gain in turbine efficiency will be 
reflected in a reduced propellant consumption, thus in an 
increase in payload. However, higher turbine performance 
usually entails multistage designs, which are heavier. A 
proprietary weight correlation is employed to estimate the 
impact of the turbine choices on the overall turbopump*' 
unit weight since a slowly rotating turbine is not only 


heavy by itself but also imposes a significant weight 
penalty on the pump side. An equation expressing the 
relationship between these opposing effects will be 
employed as a criterion to guide the optimization task. As 
will be presented in detail later, this composite objective 
function describes essentially the payload increment 
versus turbopump efficiency and weight. It is developed 
based on mission profile studies, engine balance 
perturbation and some detailed turbopump layout and 
stress information gained from other proprietary 
programs. 

2, APPROACH 

The overall approach to determine the optimum design is 
shown in Figure 1. The RSM is used to model the 
relationship between design variables and 
objective/constraint functions in function approximation 
stage of the overall approach. 

The approach of RSM is to perform a series of 
experiments, based on numerical analyses, semi- 
empirical formulas, or experimental testing, for a 
prescribed set of design points, and to construct a global 
approximation of the measured quantity over the design 
space. In this effort, numerical analysis is based on the 
aerodynamic design software. Meanline Flow Path 
Generator, for rapid analyses of turbine How' fields. 
Using overall turbine and stage input, the Meanline code 
first generates a candidate turbine flow path and displays 
a plot of the elevation view. The code, then, runs a 
meanline analysis, calculating gas conditions, velocity 
triangles, and required number of airfoils, predicted 
efficiency and power output. A calculation to predict 
turbine weight is also included. The run time required for 
analysis of a three-stage turbine on the current version of 
the code is less than one second on a Pentium II PC. The 
MUirces of performance losses due to airfoil profile, 
secondary end wall, trailing edge blockage, trailing edge 
shock, leading edge shock and unshrouded blade tip 
leakage are also included. 

Second-order polynomials are used for the response 
surface approximations for which extremal points are 
easily found by standard constrained optimization 
algorithms. The main advantages of RSM over other 
optimization tools such as gradient-based search 
algorithms are that it requires minimal interfacing with 
the analysis tools and avoids the need for expensive 
derivative calculations. To construct a second-order 
polynomial of A r design variables, the number of 
coefficients to be fixed are (N+l )t N+2)/2. 

As a first step of the overall approach of Figure i, a 
generic design box obtained by coding all design variable 
io the range i-l.+ l) is considered. Coding is based on the 
following formula: 
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- ( X max, ~ x min) ( X imx + X min . ) 

X : = X: : — + (A) 

1 1 2 

where xp real value of the design variable 

xj : normalized value of the design variable 
x^: maximum real value of the design variable 
x mio : minimum real value of the design variable 

The coding requires the information of the maximum and 
minimum values for each design variable that can be 
obtained from Tables 1-3. 


55 vv : total sum of squares about the mean 


55 V>1 = SS R + SS r (6) 

R 2 measures the proportion of the variation in the 
response around the mean that can be attributed to terms 
in the model rather than to random error 3 . R 2 a is an R 2 
value adjusted to account for the degrees of freedom in 
the model and is given by 


,= 55 g /(it z p) asi _ 
55 /(/i- 1) 


f n~' ' 


(i -/n 


n-p 


(7) 


The response surfaces of this study are generated by 
standard least-squares regression using JMP 2 , a statistical 
analysis software having a variety of statistical analyses 
functions. The global fit and prediction accuracies of the 
response surfaces are assessed through statistical 
measures such as the t-statistic, or t-ratio. nns-error, 
variation 1 1 \ The t-statistic is determined by 


vet b ) 

where hr. regression coefficient 

xe(b } y. standard error of the regression coeltieient 
and it is given by. 


Since R 2 will always increase as terms are added to the 
model, the overall assessment of the model may be better 
judged from R] . 

The polynomial-based RSM techniques are effective in 
representing the global characteristics of the design 
space. It can filter the noises associated with individual 
design data. On the other hand, depending on the order of 
polynomial employed and the shape of the actual 
response surface, the RSM can introduce a substantial 
error in certain region of the design space. An 
optimization scheme requiring large amounts of data and 
a large evaluation time to generate meaningful results is 
hardly useful. 


\c( b t — V C .. 1 3 ) 

where C,: diagonal element of i X X ) corresponding 

to b, (X is an nxp matrix ot the levels ol the 
independent variables where n is the number ot 
'hservations and p is i he number of terms in the 
model i 

s : unbiased estimator ot the standard deviation 
of the observations and unbiased estimator of the 
nns-error in prediction based on the response 
suit ace and it is uiven hv 



where e { : difference between the observation, y,-, and the 
fitted value, y. . 

The R~ \alue is determined by 
55, v 55 vv 

where 55/r : surrfof squares of the residuals or errors 
55*: sum of squares due to regression 


The optimization technique follows qualitatively that 
adopted previously for optimizing fluid machinery such 
as diffuser, injector, and airfoil, as presented in Refs. 6 ' 9 . 
The effect of numerical noise and the interaction between 
IT'D models and RSM are addressed by Madsen et al°. In 
the paper by Tucker et al. . a first effort is made to apply 
RSM for injector optimization. Pupila et alf investigated 
: he el feet of data size and relative merits between RSM 
and neural networks in handling varying data 
characteristics. The neural network technique and the 
RSM are integrated to offer enhanced optimization 
capabilities by Shyv et alf. A main focus in the present 
work is the interplay between the number of design 
variables and the predictive capability and input 
requirement of the RSM. To ascertain required predictive 
capability of the RSM, a two-level domain refinement 
approach has been adopted. As will be demonstrated in 
the following, the accuracy of the predicted optimal 
design points based on this strategy is satisfactory. 

3. DESIGN OF EXPERIMENTS (DOE) 

The response surface method is a collection of statistical 
and mathematical techniques useful for developing, 
improving, and optimizing processes and this provides an 
overall perspective of the system response within the 
design space 1 . The representation of the design space is. 
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therefore, important. In order to help to minimize the 
effect of noise on the fitted polynomial, and to improve 
the representation of the design space, design of 
experiments (DOE) procedure can be used. There are a 
number of different DOE techniques reported in the 
literature 10 17 . For example, Unal et al 10 discussed the face 
centered composite designs and D-optimal designs for 
representation of the design space for wing-body 
configuration of a launch vehicle. They showed that D- 
optimal design provides an efficient approach for 
approximating model building and multidisciplinary 
optimization. Unal et al. li studied response surface model 
building using orthogonal arrays in computer experiments 
for reusable launch vehicle and illustrated that using this 
technique minimizes design, development, test and 
evaluation cost. Similar results were obtained using three 
level fractional factorial experimental design 1 Unal and 
Dean 13 studied the robust design method based on the 
Taguchi method 14 * 1 ^ to determine the optimum 
configuration ot design parameters lor performance, 
quality and cost. They demonstrated that using such a 
robust design method for selection of design points is a 
systematic and efficient approach for determining the 
optimum configuration. 

The face centered composite design (FCCD) creates a 
design space composed of eight corners ot the cube, lour 
center of faces and the center ol the cube, figure 2 shows 
face centered composite design points tor three design 
variables. The FCCD yields (2 N + 2A/+/> points, where A 
is the number of design variables. It is more effective 
when the number of design variables is modest, say. no 
larger than 5 or 0. The FCCD is widely used for fitting 
second-order response surface 1 . 

A D-Optimal design minimizes the generalized variance. 

nf the estimates which is equivalent to maximizing 
[he determinant of the moment matrix. W . 

I I | X ' X | .QV 

I M = 18) 

l n p 

P he D-Optimal design approach requires the knowledge 
of the properties of polynomial model in selecting the 
design points. 

An orthogonal array (OA) is a fractional factorial matrix 
that assures a balanced comparison of levels of any factor 
or interaction of factors. Because the points are not 
necessarily at vertices, the orthogonal array can be more 
robust than the face centered cubic design. Based on the 
desiszn of experiments theory. OA can significantly 
reduces the number of experimental configurations. 

In this study, although the majority of the work is based 
on the FCCD approach, alternative representations of the 
design space are performed by using D-Optimal design 


and OA design. We have considered I-, 2- and 3-stage 
turbine. There are 6 design parameters for single- stage 
turbine case chosen as the mean diameter, RPM y blade 
annulus area, vane axial chord, blade axial chord, and 
stage reaction. For 2-stage turbine, mean diameter, RPM , 
exit blade annulus area, 1 st blade height (% of exit blade), 
i sl vane axial chord, 1 st blade axial chord, 2 nd vane axial 
chord, 2 nd blade axial chord, 1 st stage reaction, 2 nd stage 
reaction, and 1 st stage work fraction, are chosen and there 
are, in total, 11 design parameters. There are 15 design 
parameters for 3-stage case determined as mean diameter, 
RPM , exit blade annulus area, U c blade height (% of exit 
blade), 1 st vane axial chord. 1 st blade axial chord, 2 nd vane 
axial chord, 2 nd blade axial chord, 3 rd vane axial chord, 
3 rd blade axial chord, 1 st stage reaction, 2 nd stage reaction, 
3 rd stage reaction, I st stage work fraction, 2 nd stage work 
fraction, and 3 rd stage work fraction. Table L-3 show the 
maximum and minimum values of these parameters as 
well as their baseline values. 

With 6-input parameters of single-stage turbine, FCCD 
produces 77-data, but Meanline code produced results for 
76-data excluding one unrealistic case. Therefore, 76- 
data is used to approximate the single-stage turbine 
characteristics. With 11-input parameters of 2-stage 
turbine, FCCD yields 2,071 -data, but Meanline code 
produced results for 1990-data and this set is used to 
approximate the 2-stage turbine characteristics. For 3- 
stage having 15 -input parameters, FCCD creates 32,799- 
data based on the formula of 2' N + 2.V+/ with N-I5 
demonstrating the curse of dimensionality. For such 
cases, a statistical method can be applied to reduce the 
number of data in an efficient way. In this effort, to 
reduce the data set of 3-stageturbine of 15-dimension, D- 
Optimal design is adopted to minimize the generalized 
variance of the estimates. With the D-optimal criterion, 
the number ot data is reduced to 2500; the Meanline code 
worked for 2235 of them. In this paper, the orthogonal 
arrays are also applied for the 2-stage turbine case using a 
public domain software developed by Owen 18 . 

4. THE OPTIMIZATION PROBLEM 

The equation describing the response as given by JMP is 
input to Excel Solver. Solver is an optimization toolbox 
included with Microsoft Excel which uses the 
Generalized Reduced Gradient method 4 5 to find the 
maximum or minimum of a tunction with given 
constraints. 

The optimization problem at hand is a constrained 
optimization problem, which can be formulated as 
min{f( x )} subject to lb < x < ub. where lb is the lower 

boundary vector and ub is the upper boundary vector of 
the design variables vector x. Since the goal is to 
maximize objective tunction therefore fix) can be written 
as -gix). where g(x) is the objective function. Minimizing 


4 
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of fix) gives the same solution as maximizing the 
objective function g(x). Additional linear or nonlinear 
constraints can be incorporated if required. 

In this study, the purpose is maximizing the turbine 
efficiency, 0, and minimizing the overall weight, VV 7 , 
simultaneously. The response surface method can handle 
a multi-criteria optimization task in a straightforward 
manner by building a composite response surface from 
individual response surfaces. This composite response 
surface is referred to as the desirability function. The 
desirability function for r). i.e., tV/, can be defined which is 
to be maximized as 


5, RESULTS AND DISCUSSION 

The properties of the response surfaces obtained for rj, 
W, and Apay are shown in Tables 4 and 5 for l-,2-, and 3- 
stage turbine designs. 

For the single-stage case, there are 28-unknown 
coefficients needed for constructing the 2 nd -order 
response surface, 78 for the 2-stage and 136 for the 3- 
stage case. The quality of the fit can be evaluated by 
comparing the adjusted root mean square error (rms- 
error) shown in Table 5. 


d, = 11 T1 ™- (9) 

k "Hnux ^1 nun j 

and the desirability function for W. i.e.. d :% can be defined 
which is to be minimized as 


W - W 

:im 

w - w 

[mil ::ux 


I 10) 


where powers s and t are weighting factors which are set 
according to the role of the response in composite 
desirability function, i.e.. d. defined as follows: 


d = v'd, ' d : (11 » 

Another wav of finding optimum values of tj and \Y 
simultaneously is to maximize payload increment, Apay, 
which is a function of these two parameters in the 
following manner. 

Apay =e,xl00x< r)-r| h H W-W,j < 12) 

where r} b is the baseline efficiency and W b is the baseline 
weight. Apay function represents the amount of increase 
in payload capacity. The results of both payload 
increment based and composite desirability function 
based optimization are illustrated for 1. 2, and 3-stage 
designs. For the composite desirability function based 
optimization, different combinations of the power of d { 
and d 2 are considered with different values of t and s. 


RSM-based approximations together with the Excel 
Solver to find the maximum or minimum of the objective 
function with a given constraints is used to find the 
optimum point are obtained for all cases. Different 
starting points are tried to avoid local maximum and the 
optimum values of p, W and Apay with the corresponding 
design parameters are determined. Table 6 shows the 
optimum values of ?], W and Apay calculated for both 
Apax (Eqn. 12) based optimization and composite 
desirability function based optimization (Eqn. 11) for 
(7=/. s=0), ( t-0 , s=ll and (t=l, s=l) cases for single- 
stage turbine. The case (t=f s=0) represents the 
optimization based on weight only, whereas (7=0, s=I) 
represents the optimization based on efficiency only. The 
results shown in this table are comparable with the 
corresponding Meanline runs with the highest error of 
57c for Apay for single-stage turbine. Table 7 show's the 
optimum values of p, W and Apax calculated for both 
Apay based optimization and composite desirability 
function based optimization for (/=/. .v=0), (7=0. 5=/J, 
and (/=/, s=l) cases for the 2-stage turbine design.. 
When the results shown m this table are compared with 
the corresponding Meanltne runs of the same design 
parameters, it is observed that the percentage error is 
increased up to 13.5% for Apay. Table 8 shows the 
optimum values of rj, W and Apay calculated for the 3- 
stage turbine with the same approach. When the results 
shown in this table are compared with the corresponding 
.Meanltne runs, it is observed that the percentage error is 
increased up to 14.6 % for Apax indicating that the 
accuracy of the response surfaces constructed for this 
case is poor. 


The pitchline speed, V puch , and the lumped inertia 
measure. AN\ are used as the design constraints when 
finding the optimum solutions. 

V jlIh;ll = DxRPM I 13) 

AN 2 =AannxRPM : (14) 

where D is meanline diameter. * 


Because the accuracy of the response surface is less than 
satisfactory for 2- and 3-stage cases, we have reduced the 
size of the parameter space with the intention of 
improving the fidelity of the response surface. The details 
of the 1/5 reduced design spaces are shown in Table 9 - 
12. The new' design space is based on the optimal values 
identified for r] & W ((7=/, 5=0), (7=0, s=l) & (7=7, 

', = //) and Apay based optimization cases. With these 
refined designed spaces, substantial improvement of the 
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response surface fit accuracy is observed for 2-stage and 
also for 3-stage (Tables 13 and 14). 

Based on the results obtained, the following observations 
can be made: 

(i) To ascertain required predictive capability of the 
RSM, a two-level domain refinement strategy 
has been adopted. The accuracy of the predicted 
optimal design points based on this approach is 
shown to be satisfactory. 

(ii) For Apay-based optimization, the 2-stage turbine 
gives the best Apay result. As the number of the 
stage increases, we see that efficiency improves 
while the weight increases also. According to the 
formula for Apay, the improvement in efficiency 
can't compensate the penalty from higher 
weight. As shown in Figure 3, the mean 
diameter, speed, and the exit blade area exhibit 
distinct trends. Specifically, the diameter 
decreases, .speed increases, and annulus area 
decreases with increasing number of stages. It is 
interesting to observe that none of these design 
parameters are toward the limiting values listed 
in Tables 1-3. indicating that the optimal designs 
result from compromises between competing 
parametric trends. For such cases, a tormal 
optimizer such as the present response surface 
method is very useful. 

ml) For both weight only optimization </=/, a -0) 
and efficiency only optimization ir=U, ,v=/h as 
expected, the single-stage design gives the 
smallest weight. The efficiency also improves as 
the number of stages increases (see fables 15 
and 16. and Figure 4 for summaries). It is 
interesting to see that in both cases, the 
-elections ol the diameter, speed and annulus 
area are insensitive to the number of stages. For 
example, for weight only optimization it- 1, 
s=0). D approaches minimum, and RPM 
approaches maximum, while the annulus area is 
governed by the design constraint between it and 
RPM. 

i iv) On the other hand, for the efficiency only 
optimization it—O. \ = if . for all three stage 
designs, the annulus area approaches maximum 
and RPM is governed by the design constraint. 
The diameter, D, on the other hand, takes middle 
values within the design range (Figure 5). 
iv) Similar to the one based on A pay. the case of 
if=/. s=i ) is also a compromise between weight 
and efficiency. However, because of the different 
mathematical formulas adopted, different design 
selections than the ones based on Apay result. It 
seems clear that the precise definition of the 
optimal criterion substantially influence the 
selection of the optimal design. 


In order to determine how the optimum solution for Apay 
changes as the weighting constants, t and s , are changing, 
three designs of (7=7, s=0), ( t-0 , s-1) and (7=7, s=l) 
denoted as £2 { ,Q 2 and Q 3 respectively, are selected as 

candidate vertex points to define a plane referred as a- 
Plane as follows. 

Q = a jQj + a 2 Q. 2 + a 3 Q 3 (15) 

where = 1 and 0 < a l < 1 for /=l, 2, 3 

Equation 12 was used to obtain 66 design points 
distributed on the alpha plane for 2-stage turbine. Figure 
6 illustrates the contour plots over the a-Plane for 
optimum Apay. This figure shows that Apay reaches its 
maximum value for a range of t. 0<t<I for s=J . 

6. ORTHOGONAL ARRAYS FOR 2-STAGE 
TURBINE DESIGN 

Although the majority of the present work is based on the 
FCCD approach, orthogonal arrays are constructed to 
investigate the efficiency of orthogonal array designs in 
representing the design space for 2-stage turbine. For this 
purpose. 249 design points are selected using OA 
designs, fable 17 shows the comparison of the statistics 
of the second-order response surfaces generated for r], \V 
and Apay by using 1990-data generated by face centered 
composite design and 249-data selected by OA method. 
This table illustrates that the fidelity of the response 
surface generated for design space of 249 data, based on 
orthogonal arrays, are comparable with that of 1990 data 
based on the face centered criterion. The response surface 
models are also assessed by using 78-test data to 
determine the predictive accuracy of these models. Table 
IS presents that the testing nns-crrors of response 
surfaces generated are 1.65% for r) and 0.96% for W 
using 249-data, and 1.67% for rf and 1,21% for Housing 
1990-data. The results of optimization based on Apay and 
composite desirability function of 7] & W with 249-data 
selected by orthogonal arrays are shown in Table 19 for 
t he original design space and in Table 20 for the refined 
design space. When these results are compared with the 
results of 1990-data presented in Tables 7 and 13, it is 
observed that the optimum rj, W and Apay are largely 
consistent. However, it is also observed from Figure 7 
which shows the comparison of the design variables for 
optimization based on {Apay). some of the design 
variables are different even though optimum rj, W and 
Apay are consistent. This shows that there are multiple 
points in the design space which yield comparable 
performance. Nevertheless, it remains true that the two- 
- stage turbine is most suitable from a payload point of 
view. 


6 
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7. CONCLUDING REMARKS 

In summary, the result of the RSM computation indicates 
that indeed the efficiency rises quickly from 1 stage to 2 
stages but the increase is much less pronounced with 3 
stages. A 1 -stage turbine performs poorly under the 
engine balance boundary condition. A significant portion 
of kinetic energy is lost at the turbine discharge of the 1- 
stage design due to high pressure ratio and high energy 
content, mostly hydrogen, of the working fluid. Adding a 
2 nd stage recovers most of that wasted energy resulting in 
much better efficiency for a 2-stage turbine. An extra 3 fd 
stage only improves the efficiency slightly. 
Understandably, the turbopump weight also increases 
substantially from 1 to 3 stages even though the 3-stage 
turbine diameter is smaller. The smaller diameter is the 
direct outcome of higher RPM that is the result of lower 
exit annulus area satisfying the constraint defined tor AN". 
The exit annulus area is smaller because of the lower 
pressure ratio per stage. However, the 3-stage turbine is 
much longer and requires a change in bearing 
configuration that adds significantly into the overall 
weight. The optimum 2-stage turbine resulting from the 
RSM optimization is consistent with a design produced by 
an experienced engineer; namely, that most of the work is 
done by the U‘ stage at very low reaction. By varying 
from 1- to 2-to 3-stages, we observe that while the size of 
the training data increase naturally with the number ot 
design variables, the actual need is case dependent, 
furthermore, it seems that the selection ot the data 
distribution can be more critical than the data size. Present 
investigation has also demonstrated that the criteria tor 
selecting the data base exhibit significant impact on the 
efficiency and effectiveness ot the construction ot the 
response surface. 
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Table 1. Design Space for Single-Stage Turbine (All geometric design variables are normalized by the baseline 


values) 


Variable 

Lower Limit 

Upper Limit 

Mean Diameter, D 

0.50 

1.50 

Speed, RPM 

0.70 

1.30 

Blade Annulus Area. A aIin 

0.70 

1.30 

Vane Axial Chord, c, ( 

0.39 

1.71 

Blade Axial Chord, e b 

0.26 

1.14 

Stage Reaction, sr 

0.0% 

50% 


Fable 2, Design Space for 2-Stage Turbine (All geometric design variables are normalized by the baseline values) 

Variable Lower Limit Upper Limit 

1.50 
1.30 
1.30 
1 .50 
1.71 
1.14 
1.41 
1.13 
50% 

50% 

85% 


Mean Diameter, D 

0.50 

Speed, RPM 

0.70 

Blade Annulus Area. A ;mn 

0.70 

1 i! Blade Height ( % of Exit Blade), h| 

0.90 

1 st Vane Axial Chord. c vl 

0.39 

T' Blade Axial Chord. c bl 

0.26 

2 lld Vane Axial Chord, c v2 

0.21 

2 nd Blade Axial Chord, c b2 

0.17 

1 st Stage Reaction, sr t 

0.0% 

2 nd Stage Reaction, sr 2 

0.0% 

1 st Work Fraction, wf| 

50% 
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Table 3. Design Space for 3-Stage Turbine (All geometric design variables are normalized by the baseline values) 


Variable 

Lower Limit 

Upper Limit 

Mean Diameter, D 

0.50 

1.50 

Speed, RPM 

0.70 

1.30 

Blade Annulus Area, A^ n 

0.70 

1.30 

I st Blade Height (% of Exit Blade) , hj 

0.90 

1.50 

1 st Vane Axial Chord, c vl 

0.39 

1.71 

1 st Blade Axial Chord, c M 

0.26 

1.14 

2 nd Vane Axial Chord. c v2 

0.21 

1.41 

2 nd Blade Axial Chord, c b2 

0.17 

1.13 

3 rd Vane Axial Chord, c v3 

0.21 

1.41 

3 rd Blade Axial Chord, cv, 

0.17 

1.13 

1 st Stage Reaction, sf| 

0.0% 

50% 

2 nd Stage Reaction, sr 2 

0.0% 

50% 

3 rd Stage Reaction, sr 3 

0.0% 

50% 

1 st Work Fraction, wf } 

40% 

80% 

2 nd Work Fraction, wf 2 

30% 

10% 


Table 4. Response Surface Summary for 1, 2 and 3-Stage Turbine 



No. of Design Parameters 

No. of Coefficients 

No. of Design Points 

1 -Stage 

6 

28 

76 

2-Stage 

1 1 

78 

1990 

3-Stage 

15 

136 

2235 


Table 5. The quality of the Second-Order Response Surface obtained for rj, W, and Apay 1 , 2 and 3-Stage Turbine 
(Mean values of )). IV and Apay are normalized by the baseline values) 




0 

W 

Apay 


R z 

0.998 

0.999 

0.998 


Ra : 

0.997 

0.999 

0.997 

1 -Stage 

mis- error 

2.509? 

0.82% 

4.09% 


Mean 

0.57 

0.60 

-0.43 


R : 

0.995 

0.996 

0.995 


Ra : 

0.994 

0.996 

0.995 

2-Stage 

mis- error 

1.31% 

2.56% 

9.58% 


Mean 

0.78 

0.86 

-0.24 


R‘ 

0.989 

0.989 

0.994 


Ra : 

0.988 

0.988 

0.994 

3 -Stage 

mis- error 

2.0% 

2.9% 

8.4% 


Mean 

0.89 

1.41 

-0.26 


9 

American Institute of Aeronautics and Astronautics 




AIAA-2000-3242 


Table 6. Optimization based on Apay and composite desirability function of T] and W for single-stage turbine 
(All geometric design variables and output parameters are normalized by the baseline values) 




Wop, 

Apay opl 

D 

RPM 

A-ann 

Cy 

c b 

sr 

RSM(Zlpav) 

0.766 

0.731 

-0.214 


0.975 

1.166 

1.706 

0.880 

0 

Meanline 

0.797 

0.733 

-0.193 

Error % of mean 

2.9 

0.3 

CO 

A 






RSM <7=7, s=0) 

0.399 

0.407 

-0.611 

0.502 

1.284 

0.699 

0.394 

0.264 

0.5 

Meanline 

0.383 

0.402 

-0.623 

Error % of mean 

1.8 

l 

2.7 






RSM(7=0. s=l) 

0.781 

0.762 

-0.216 

1.260 

0.915 

1.300 

1.575 

0.880 

0 

Meanline 

0.797 

0.762 

-0.199 

Error % of mean 

2.2 

0.03 

3.8 






RSM(t=I, s=J) 

0.702 

0.583 

-0.261 

0.895 

1.284 

0.699 

1.706 

0.264 

0 

Meanline 

0.718 

0.588 

-0.239 

Error % of mean 

3.1 

0.8 

s i 



Table 7. Optimization based on Apay and composite desirability function ot t] and \V tor 2-stage turbine tor 
original design space (All geometric design variables and output parameters are normalized by the baseline values) 



Oopt 

W 0pt 

Apayop, 

D 

RPM 

•^ann 

h, 

C v i 

C v 7 

c bl 

Cb2 

sr i 

sr 2 

Wfi 

RSM (Apay) 

1.10 

1.05 

0.11 

1.16 

0.99 

1.14 

1.50 

1.57 

0.97 

0.71 

0.68 

0 

0.50 

0.85 

Meanline 

1.12. 

1.05 


Error kof mean 

IfeMil 








RSM(7=/. v =0) 

0.65 



0.50 i 1.27 

L 

0.70 

1.50 

1.71 

1.76 

0.92 

1.13 

0.50 

0.50 

0.50 

Meanline 

0.65 

0.65 


Error 9k) f mean 

1.10 









— 


RSM U=(). s-i ) 

1.10 



1.24 ; 0.92 

________ ! i 

1.30 

1.50 

1.44 

1.06 

0.71 

0.79 

0 

0 

0.80 



Error 7cof mean 





1 r r 



RSM<r=7. v = / > 

0.99 



0.91 | 1.27 | 0.70 | 1.50 j 

1.71 

1.76 j 

0.85 j 

1.13 

0 

0 j 0.50 

Meanline 0.99 - 0.84 

0.03 

Error kof mean 

0.30 

1.70 

0.40 



Table 8. Optimization based on Apax and composite desirability function of i] and W for o-stage turbine for 
original design space (All geometric design variables and output parameters are normalized by the baseline values) 




W 

vv opt 

Apay opl 


RSM (Apay) 

1.24 

1.62 

0.14 


Meanline 

1.21 

1.57 

0.10 


Error %of mean 

4.72 

3.52 

14.59 


RSM (t=l, s=0) 

0.85 

1.13 

-0.23 


Meanline 

0.83 

1.13 

-0.26 


Error k of mean 

2.13 

0.34 

14.39 


RSM (t=0. s=!) 

1.26 

1.74 

0.12 


Meanline 

1.23 

1.69 

0.09 


Error %of mean 

4.27 

3.10 

13.54 


RSM (t=L s=I) 

1.08 

1.33 

0.05 


Meanline , 

1.10 

L34 

0.04 


Error 7c of medn 

2.16 

0.59 

4.11 
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D 

RPM 

Aann 

hi 

C v i 

Cy2 

Cv3 

Cbl 

C b 2 

c b3 

sri 

sr 2 

srj 

wf. 

wf 2 

RSM ( Apa\ ) 

1.07 

1.07 

0.98 

1.50 

1.59 

1.09 


0.56 

1.02 

0.89 

0 

0.5 

0.5 

0.6 

0.3 

RSM (t=l, s=0) 

0.50 

1.28 

0.70 

0.90 

1.71 

1.76 

1.41 

0.73 

1.41 

0.17 

0.5 

0.5 

0.5 

0.4 

0.3 

RSM (t=0, s=l) 

1.19 

0.96 

1.20 

1.50 

1.44 

1.06 

0.78 

0.56 

0.99 

0.90 

0 

0.5 

0.5 

0.6 

0.3 

RSM (t=l, s=l ) 

0.91 

1.29 

0.70 

1.50 

1.71 

0.62 

0.78 

0.17 

0.21 

1.13 

0 

0 

0 

0.6 

0.1 


Table 9. Upper and Lower Limits of the Design Parameters of the refined designed spaces for 2-stage turbine 
(All geometric design variables are normalized by the baseline values) 




D 

RPM 

A 

^ann 

hi 

C v i 

C v 2 

c bl 

c b2 

sri 

sr 2 

Wfi 


Max 

1.50 

1.30 

1.30 

1.50 

1.71 

1.76 

0.92 

1.13 

0.50 

0.50 

0.85 

Original 

Min 

0.50 

0.70 

0.70 

0.90 

0.39 

0.26 

0.21 

0.17 

0.00 

0.00 

0.50 

Refined 

Max 

1.26 

1.05 

1.20 

1.50 

1.71 

1.11 

0.79 

0.79 

0.05 

0.50 

0.85 

(Apay) 

Min 

1.06 

0.93 

1.08 

1.44 

1.44 

0.81 

0.64 

0.60 

0.00 

0.45 

0.82 

Refined 

Max 

0.60 

1.30 

0.76 

1.50 

1.71 

1.76 

0.92 

1.13 

0.50 

0.50 

0.54 

(t-I, s=0) 

Min 

0.50 

1.21 

0.70 

1.44 

1.57 

1.61 

0.85 

1.03 

0.45 

0.45 

0.50 

Refined 

Max 

1.34 

0.98 

1.30 

1.50 

1.63 

1.20 

0.76 

0.90 

0.05 

0.05 

0.82 

(t=0, 5= J ) 

Min 

1.14 

0.86 

1.24 

1.44 

1.36 

0.91 

0.62 

0.70 

0.00 

0.00 

0.75 

Refined 

Max 

1.00 

1.30 

0.76 

1.50 

1.71 

1.76 

0.92 

1.13 

0.05 

0.05 

0.54 

(t-i , s= I ) 

Min 

0.80 

1.21 

0.70 

1.44 

1.57 

1.61 

0.80 

1.03 

0.00 

0.00 

0.50 



Table 10, Center Coordinates of the refined designed spaces for 2-stage turbine 





(All geometric design 

variables are normalized by the baseline values) 





D 

RPM 

Amn 

h ( 

C V 1 

Cv2 

Cbl 

Cb2 

srj 

sr 2 

Wfi 

Original 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Refined {Apay) 

1.16 

0.99 

1.14 

1.50 

1.57 

0.97 

0.71 

0.68 

0 

0.5 

0.85 

Refined (t=l. 

s=0) 

0.50 

1.27 

0.70 

1.50 

1.71 

1.76 

0.92 

1.13 

0.5 

0.5 

0.5 

Refined (t=(), 

s=l) 

1.24 

0.92 

1.30 

1.50 

1.50 

1.06 

0.69 

0.80 

0.0 

0.0 

0.78 

Refined (t=I, 

s=n 

0.91 

1.27 

0.70 

1.50 

1.71 

1.76 

0.85 

1.13 

0 

0 

0.50 


Table 11. Upper and Lower Limits of the Design Parameters of the refined designed spaces for 3-stage turbine 
(All geometric design variables are normalized by the baseline values) 




D 

RPM 

A 

^ann 

h, 

C V [ 

C v2 

C V 3 

WEM 


B39 


EE 

ESI 

ESI 

ii 


Max 

1.50 

1.30 

1.30 

1.50 

1.71 

1.76 

1.41 

0.73 

1.41 

1.13 

0.5 

0.5 

0.5 

0.8 

0.3 

Original 

Min 

0.50 

0.70 

0.70 

0.90 

0.39 

0.26 

0.21 

0.17 

0.21 

0.17 

0.0 

0.0 

0.0 

0.4 

0.1 

Refined 

Max 

1.17 

1.13 

1.05 

1.50 

1.71 

1.24 

0.99 

0.62 

1.14 

0.99 

0.05 

0.50 

0.50 

0.62 

0.30 

(Apav) 

Min 

0.97 

1.01 

0.92 

1.44 

1.46 

0.94 

0.75 

0.51 

0.90 

0.79 

0.00 

0.45 

0.45 

0.56 

0.26 

Refi ned 

Max 

0.60 

1.30 

0.76 

0.96 

1.71 

1.76 

1.41 

0.73 

1.41 

0.26 

0.50 

0.50 

0.50 

0.44 

0.30 

(t-I , s-0) 

Min 

0.50 

[ ii 

0.70 

0.90 

1.57 

1.61 

1.29 

0.68 

1.29 

0.17 

0.45 

0.45 

0.45 

0.40 

0.28 

Refined 

Max 

1.29 

1.02 

1.26 

1.50 

1.60 

1.18 

0.93 

0.60 

1.10 

0.99 

0.05 

0.50 

0.50 

0.62 

0.30 

n 

ii 

Min 

1.09 

0.90 

1.14 

1.44 

1.34 

0.88 

0.69 

0.49 

0.86 

0.80 

0.00 

0.45 

0.45 

0.56 

0.26 

Refined 

Max 

1.00 

1.30 

0.76 

1.50 

1.71 

0.74 

0.88 

0.23 

0.33 

1.13 

0.05 

0.05 

0.05 

0.64 

0.14 

a-L s-I ) 

Min 

0.80 

1.23 

0.70 

1.44 

1.57 

0.44 

0.64 

0.17 

0.21 

1.03 

0.00 

0.00 

0.00 

0.58 

0.10 
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Table 12. Center Coordinates of the refined designed spaces for 3-stage turbine 
(All geometric design variables are normalized by the baseline values) 



D 

RPM 

Ajuuj 

h, 

Cv| 

C v2 

C v 3 




sr i 

sr 2 

Sf 3 

wfi 

wf 2 

Original 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Refined (Apay) 

1.07 

1.07 

0.98 

1.50 

1.57 

1.06 

0.85 

0.56 

0.99 

0.90 

0.0 

0.5 

0.5 

0.6 

0.3 

Refined (t-I t s=0) 

0.50 

1.28 

0.70 

0.90 

1.71 

1.76 

1.41 

0.73 

1.41 

0.17 

0.5 

0.5 

0.5 

0.4 

0.3 

Refined (t=0, s=J) 

1.19 

0.96 

1.20 

1.50 

1.44 

1.06 

0.78 

0.56 

0.99 

0.90 

0.0 

0.5 

0.5 

0.6 

0.3 

Refined (7=/, s-I ) 

1.00 

1.00 

1.00 

1.20 

1.05 

1.06 

0.85 

0.45 

0.85 

0.68 

0.3 

0.3 

0.3 

0.6 

0.2 


Table 13. Optimization based on Apay and composite desirability function of rj and W for 2-stage turbine for 
refined design space (All geometric design variables and output parameters are normalized by the baseline values) 



TW 

W* 

Apayopt 

D 

RPM 

Ajpm 

ID 


C v 2 

Cbl 

c b2 

sri 

sr 2 

Wfi 

RSM (Apay) 

1.13 

1.04 

0.15 

1.12 

1.02 

1.08 

1.5C 



0.71 

0.62 

0.1 

0.5 

0.9 

! Meanline 

1.13 

1.04 

0.15 

Error 7c of mean 

0.03 

0.02 

0.16 






■Hi 



RSM (t=L s=0) 

0.65 

0.65 

-0.35 


1.27 

0.70 

1.50 

1.71 

1.76 

0.92 

j 1.13 




Meaniine 

0.65 

0.65 

-0.35 

Error %of mean 

0.00 

0.00 

0.01 



RSM(t=0, s=l) 

1.15 

1.10 

0.15 

1.23 


1.30 

1.50 

1.31 

0.88 

0.71 

0.73 

0.1 

0 

0.8 

Meanline 

1.15 

1.10 

0.15 

Error 7coi mean 

0.02 

0.01 

0.09 






RSM(/=/, .v=/) 

1.00 

0.85 

0.04 


1.27 

0.70 

1.50 

1.71 

1.58 

0.92 j 

1.01 

0 

0 

0.5 



0.04 






Table 14. Optimization based on Apay and composite desirability function ot rj and If tor 3-stage turbine tor 
refined design space (All geometric design variables and output parameters are normalized by the baseline values) 



"Hop! 

W 

yy opl 

Apay op , 

RSM ( Apav ) 

1.20 

1.54 

0.11 

Meanline 

1.21 

1.54 

0 . 1 1 

Error %of mean 

0.22 

0.13 

0.41 

RSM (t=l, s=0) 

0.82 

1.13 

-0.27 

Meanline 

0.82 

1.13 

-0.27 

Error 7caf mean 

0.04 

0.02 

0.31 

RSM (t=0, s=l) 

1.24 

1.75 

0.09 

Meanline 

1.23 

1.72 

0.09 

Error %of mean 

1.39 

1.29 

1.35 



D 

MM 


h| 

. .Cirlifei 


Cy3 


Bil 

EH 

sri 


E?1 


MM 

RSM (Apay) 

1.03 

1.11 

0.92 

1.50 

1.46 

0.94 

0.75 

0.51 

0.90 

0.79 

0.04 

0.5 

0.5 

0.62 

0.28 

RSM(r=/, s=0) 

0.50 

1.27 

0.70 

0.96 

1.71 

1.76 

1.41 

0.73 

1.41 

0.17 

0.5 

0.5 

0.5 

0.4 

0.3 

RSM (t=a s=I ) 

1.20 

0.94 

1.26 

i .47 

1.35 

0.88 

0.69 

0.54 

0.86 

0.80 

0 

0.45 

0.45 

0.6 

0.3 
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Table 15. Optimization summary for 1. 2 and 3-stage turbine with response surface in original design space (All 

output parameters are normalized by the baseline values) 




* 1 opt 

Wop, 

Apay,*, 


1 -stage 

0.77 

0.73 

- 0.21 

Apay 

2-stage 

1.10 

1.05 

0.11 


3 -stage 

1.24 

1.62 

0.14 


1 -stage 

0.40 

0.41 

-0.61 

(t=l, s=0) 

2-stage 

0.65 

0.66 

-0.34 


3-stage 

0.85 

1.13 

-0.23 


1 -stage 

0.78 

0.76 

-0.22 

(t=0, s=l) 

2-stage 

1.10 

1.10 

0.11 


3-stage 

1.26 

1.74 

0.12 


l -stage 

0.70 

0.58 

-0.26 

(t=l. s=l) 

2-stage 

0.99 

0.85 

0.03 


3-stage 

1.08 

1.33 

0.05 


Table 16. Optimization summary for 1. 2 and 3-stage turbine with response surface in refined design space 
(All output parameters are normalized by the baseline values) 




^lopt 

Wop, 

Apayop, 


1 -stage 

0.77 

0.73 

-0.21 

Apay 

2-stage 

1.13 

1.04 

0.15 


3 -stage 

1.20 

1.54 

0.11 


1 -stage 

0.40 

0.41 

-0.61 

it=i. j ■=<)) 

2-stage 

0.65 

0.65 

-0.35 


3 -stage 

0.82 

1.13 

-0.27 


l -stage 

0.78 

0.76 

-0.22 

(!~0, S=I ) 

2-stage 

1.15 

1.10 

0.15 


3-stage 

1.24 

1.75 

0.09 


1 -stage 

0.70 

0.58 

-0.26 

u=i. s=n 

2-stage 

1.00 

0.85 

0.04 


Table 17.The quality of the Second-Order Response Surface obtained for rj. Wand A pay of 2-Stage Turbine for 

1990-data (FCCD criterion) and 249-data (OA criterion) 

(Mean values of T}, W and Apay are normalized by the baseline values) 




n 

W 

A P a y 


R : 

0.995 

0.996 

0.995 


Ra : 

0.994 

0.996 

0.995 

1990-data 

mis - error 

1.31% 

2.56% 

9.58% 


Mean 

0.78 

0.86 

-0.24 


R 2 

0.995 

0.998 

0.994 


Ra 2 

0.992 

0.997 

0.992 

249-data 

mis- error 

2.128% 

0.826%) 

20.68% 


Mean 

0.89 

0.92 

-0.11 
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Table 18. Testing of the Second-Order Response Surface obtained for r\ and W of 2-Stage Turbine for 1990-data 
(FCCD criterion) and 249-data (OA criterion) with 78- test data 


# of design points 

# of test data 

rms -error for r| (%) 

rms- error for W(%) 

249 

78 

1.65 

0.96 

1990 

78 

1.67 

1.21 


Table 19. Optimization based on Apay and composite desirability function of r| and W for 2-stage turbine for 
original design space with 249-data (OA criterion). (All geometric design variables and output parameters are 

normalized by the baseline values) 



*)oDt 

W 0pt 

Apay opt 

D 

RPM 

A 

* 'ann 

h, 

Cy 1 

C V 2 

c bl 

c b2 

sri 

sr 2 

w n 

RSM (Apay) 

1.13 



1.13 

1.02 

1.10 

1.50 

1.71 

0.44 

0.92 


0.0 

0.4 

0.9 

Meanline 

1.12 

1.03 

0.14 

Error %of mean 

1.73 

1.39 

5.90 







RSMlf=7, s=0) 

0.64 

0.64 

-0.36 

0.50 

1.27 

0.70 

0.90 

0.39 

m 


1.13 

0.5 

0.5 

0.5 

Meanline 

0.62 

0.64 

-0.38 

Error Vcof mean 

1.85 

0.62 

7.15 







RSM (t=0, s~ I ) 

1.13 

1.09 

0.14 

1.21 

0.95 

1.27 

1.50 

1.71 

0.53 

0.85 

0.68 

0.0 

0.3 

0.9 

Meanline 

1.13 

1.09 

0.14 

Error %of mean 

0.90 

1.02 

2.76 







RSM(f=/. s=n 

0.99 

0.82 

0.03 

0.91 

1.27 

0.70 

0.90 

0.39 

1.76 

0.21 

1.13 

0.5 

0.0 

0.9 

Meanline 

0.96 

0.83 

0.00 

Error v k of mean 

2.67 

0.49 

11.88 



Table 20. Optimization based on Apay and composite desirability function of r\ and W tor 2-stage turbine tor 
refined design space for 249-data (OA criterion). (All geometric design variables and output parameters are 

normalized bv the baseline values) 



Hopi 

W 

vv opt 

Apay opt 

D 

RPM 


IB 




c b2 

srj 

1 sr 2 

Wfi 

RSM {Apay) 

1.13 



1.10 

1.05 

1.03 

1.5C 

IB 



0.51 




Meanline 1.12 

1.02 

0.15 



0.01 

0.37 


mm 

HH 

1 H 

RSM(t=I. s=0) 

0.62 

0.64 

-0.38 

0.50 

1.27 

0.70 

0.90 

0.39 

0.26 

0.92 





Meanline 0.63 

0.64 

-0.38 

Error ( \ of mean 

0.0031 

0.0002 

0.0002 


IS 



si 

RSMf7=0, s=n 

1.13 

1.10 

0.15 


0.93 

1.30 

1.50 

1.57 

0.70 

0.78 




■ 

Meanline 

1.14 

1.10 

0.15 

Error %of mean 

0.005 

0.001 

0.029 




RSM|7=/, s=l) 

0.97 

0.84 

0.02 

0.90 

1.27 

0.70 

0.90 

0.39 

1.76 

0.28 

1.13 

0.5 

0.0 

0.8 

Meanline 

0.98 

0.84 

0.02 

Error r fof mean 

0.08 

0.05 
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Figure 1 . Function Approximation and Optimization Flow Chart 



Figure 2. Face Centered Composite Designs (FCCD) tor 3 Design Variables 
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Bfect of # of Turbine Stages for Optimum Bficiency 
( normalized by the baseline value) 



Optimum Diameter 



1.7 i 

i 

1.5 I 

i.3 i 
H 

1 

5 

0.9 

0.7 


Bfect of # of Turbine Stagesfor Optimum Weight 
( normalized by the baseline value) 


0.5 J- 
0 




1.20 ; 
1.15 I 
110 


cl 1.05 
QC 


1.00 
0.95 ! 
0.90 J- 


Optimum RFM 

{ normalized by the baseline value) 


* 




$ of Stage 


2 

# of Stage 


Bfect of # of Turbine Stages for Optimum Incremental 
Payload Change ( normalized by the baseline value) 


'§ 0.1 

ca 

.c 

O 


•5 0.0 I 

o 

>• 6 

CL 

75-0.1 ; 

c 

o 

2 - 0.2 i 

y 


-0.3 





# of Stage 



1.20 , 
t.15 j 
1.10 
1.05 

| 00 < 

0.95 | 

0 90 j 

0.85 I 

0.80 L 
0 


Optimum Annulus Area 
( normalized by the Daseime value) 



1 2 3 

# of Stage 


Figure 3. L: fleet of the number of turbine stage on optimum design parameters; D, RPM . and A, mn and optimum 
output parameters; 1], \V, and Apay calculated for Apay - based optimization 
(All geometric design variables and output parameters are normalized by the baseline values) 
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Bfect of # of Turbine Stages for Optirrum Bficiency 
{ normalized by the baseline value) 



# of Stage 


Bfect of # of Turbine Stagesfor Optimum Weight 
( normalized by the baseline value) 


' 1,2 T 

1.1 s 

1.0 

0.9 

■S ) 08 

| 0.7 
0.6 
0.5 

0.4 ; 

0.3 - 


<r 




0 


0.600 


0.550 


Q 0.500 


0.450 


0.400 


CL 

cr 


1.30 | 

1.29 j 

I 

1 28 
1 27 
1.26 


25 


Optimum Diameter 
( normalized by the baseline value) 



2 

# of Stage 


Optirrum RPM 

; normalized by the baseline value) 


'V. 


# of Stage 


2 

f of Stage 


4 1 


Bfect of # of Turbine Stages for Optimum Pay dad 
{ normalized by the baseline value) 


0.0 
- 0.1 
- 0.2 
■o -0.3 

O j 

S' -0.4 I 

Q_ 

-0.5 | 
- 0.6 | 
-0.7 i- 




# of Stage 


0.71 


0.70 


0.69 


Optimum Annulus Area 
( normalized by the baseline value) 


2 

# of Stage 


Figure 4. Effect of the number of turbine stage on optimum design parameters; D , RPM, and A tinn and optimum 

output parameters; rf, W, and Apay calculated 
for Weight-based optimization (/=/. s-0) 

(All geometric design variables and output parameters are normalized by the baseline values) 
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Optirrum Diameter 
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16 I 

14 | 
| 1-2 
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Bfect of # of Turbine Stagesfor Optimum Weight 
( normalized by the baseline value) 
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^ 0.05 

CL 

■ 0,10 

- 0.15 

- 0.20 
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( normalized by the baseline value) 
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1 35 
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Figure 5. Effect of the number of turbine stage on optimum design parameters; D, RPM , and A, mn and optimum 

output parameters; r\, IV, and A pay calculated 
for ?7-based optimization (t=(). s=/j 

(All geometric design variables and output parameters are normalized by the baseline values) 
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a-plot for Apay Response Surface for 2-Stage Turbine 
(normalized by the baseline value ) 



(t=l,s=0) 


(t=l.s=l) 




Dpay 
ym 0.118547 
“ 0.108076 
0.100127 
0.0825071 
0.0646329 
0.0441748 
i| 0.00650685 
^ -0.0415149 
I -0.0755047 
-0.112671 
< -0.167669 
JJ -0.224047 
B -0 279187 
™ -0.3322 


(t=i,s=0) 


Figure 6. a-Plot for Apay (normalized by the baseline values ) 
Response Surface for 2-Stage Turbine based on composite desirability 
function optimization (Effect of values of t and s on optimum Apay ) 
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(a) Original Design Space 



ib) Refined Design Space 


Figure 7. Comparison of the Design Variables for Optimization based on Payload Increment { Apay) using 1990- 
data (FCCD) and 249-data (OA) for both Original Design Space and Refined Design Space 

(DV#1: D, DV#2: RPM , DV#3: A ann , DV#4: h,, DV#5: c vl , DV#6: c v2 , DV#7: c bI , DV#8: c b2 , DV#9: sr,, DV#10: 

sr 2 , and DV#1 1: w n ) 
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